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ABSTRACT 

The cosmic abundance of the 3 He isotope has important implications for many 
fields of astrophysics. We are using the 8.665 GHz hyperfme transition of 3 He + to 
determine the 3 He/H abundance in Milky Way H II regions and planetary nebu- 
lae. This is one in a series of papers in which we discuss issues involved in deriving 
accurate 3 He/H abundance ratios from the available measurements. Here we de- 
scribe the ionization correction we use to convert the 3 He + /H + abundance, y£, to 
the 3 He/H abundance, y^. In principle the nebular ionization structure can sig- 
nificantly influence the 2/3 derived for individual sources. We find that in general 
there is insufficient information available to make a detailed ionization correction. 
Here we make a simple correction and assess its validity. The correction is based 
on radio recombination line measurements of H + and 4 He + , together with simple 
core-halo source models. We use these models to establish criteria that allow us 
to identify sources that can be accurately corrected for ionization and those that 
cannot. We argue that this effect cannot be very large for most of the sources in 
our observational sample. For a wide range of models of nebular ionization struc- 
ture we find that the ionization correction factor varies from 1 to 1.8. Although 
larger corrections are possible, there would have to be a conspiracy between the 
density and ionization structure for us to underestimate the ionization correction 
by a substantial amount. 
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INTRODUCTION 



1.1. The 3 He Experiment 

The abundance of the light isotope of helium, 3 He, is astrophysically important. Know- 
ing the 3 He/H abundance ratio can be used to test the theory of stellar nucleosynthesis; 
it gives important information needed to evaluate models of Galactic chemical evolution; 
it can help constrain Big Bang Nucleosynthesis. For over two decades we have used the 
8.665 GHz hyperfine transition of 3 He + to derive the 3 He/H abundance in the interstellar 
medium (ISM) of the Milky Way. Our 3 He sources, planetary nebulae and H II regions, 
are distributed throughout the Milky Way's disk, from the Galactic Center to the outer- 
most regions. There is no other 3 He spectral transition available that can be used to probe 
transgalactic paths with the sensitivity and accuracy of the 3 He + hyperfine line. 

We observed 3 He + using the Max-Planck-Institut fur Radioastronomie (MPIfR) 100 m 
telescope and the National Radio Astronomy Observatory (NRAO) 1 140 Foot telescope and 
Very Large Array ( VL A) . With the 100 m and VLA we primarily studied plane tary nebulae 



(Bal ser et al.lll997l . 120061 ) whereas the 140 Foot was used mostly for H II regions (IBania et al. 
19971 . hereafter Paper I). In addition to the 3 He + observations, measurements were also 
made of a variety of radio recombination line transitions (RRLs) of H, 4 He, and C, as well 
as of the thermal continuum emission from the nebulae. Paper I discussed the observations 
and the measurement errors. It is a summary of the status of the observations as of 1996 
March. We continued the 3 He + program on the 140 Foot until it was decommissioned in 
1999 July. (The last scientific spectrum observed at the 140 Foot telescope w as the 3 He 



scan that finished at 08:12:20 EDT 19 July 1999.) The final paper in this series ( IRood et al. 
20071 . Paper IV) summarizes the results of the 3 He + experiment for all the sources observed 



with the NRAO 140 Foot telescope during the period 1982-1999. Paper IV compiles the 
observed properties of the 3 He + emission and gives the final 3 He/H abundances derived for 
the NRAO 140 Foot sample of 3 He + H II regions. 

The quantity of astrophysical interest is the 3 He/H abundance ratio by number w hich 
we define as y^. This ratio provides information about both stellar (ICharbonnell Il998l ) and 
Galactic chemical evolution ( Romano et al. 2003h and con s trains cosmq l ogica l models during 
the era of primordial nucleosynthesis (jYang et al.l Il984l ; IBania et al.l 120021 ). The species 
directly accessible to observation, however, are 3 He + and H + . The collisionally excited 3 He + 
hyperfine transition directly measures the total column density of 3 He + atoms along the 
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line-of-sight within the telescope's beam. To determine the column density of H + , either 
radio continuum or hydrogen RRLs can be used. At 8.7 GHz the radio continuum emission 
in H II regions and planetary nebulae is primarily due to thermal free-free emission which 
is proportional to the emission measure, EM = J n 2 e di, where n e is the electron density 
and dt is the differential path length through the ionized nebula. The H RRL emission is 
also proportional to the emission measure. Thus neither the RRL nor the radio continuum 
data for H + are probing the total proton column density (or mass) in a straightforward 
manner. In order to determine the 3 He + /H + abundance ratio, it is therefore necessary to 
model the nebula's density structure. Models are required because in most cases detailed 
information on nebular density structure is not observationally accessible. Ne bular density 



model s and the 3 He + /H + abundance ratios derived from them are discussed by lBalser et al. 



( 119991 . hereafter Paper II). In particular, Table 5 of Paper II lists the adopted 3 He + /H + 
abundance ratios for the 21 H II regions discussed in Paper I. 

Using the observables, 3 He + and H + , to derive a 3 He/H abundance ratio, y 3 , is, however, 
a two step process. The first step is to determine the source density structure and calculate 
3 He + /H + (Paper II). The next step, the topic of this paper, is to determine the source 
ionization structure and to use this information to convert 3 He + /H + into 3 He/H. This 
requires an understanding of the ionization properties of each ionized nebula. Specifically, 
one needs to know the ionization structure of both the He and H throughout each nebula. 
Because the first ionization potentials of H (13.6 eV) and He (24.6 eV) are nearly a factor of 
two different, in principle ionization can significantly influence y% for individual sources. We 
find that in general there is insufficient information available to make a detailed ionization 
correction. Here we make a simple correction and assess its validity. 



1.2. Helium Ionization Structure in HII Regions 



Several diagnostics are commonly used to study the helium ionization structure in H II 
regions. The most direct is to measure the fraction of 4 He°, 4 He + , and 4 He ++ within the H II 
(H + ) region. Since both 4 He and 3 He have essentially the same ionization potential the 3 He + 
and 4 He + emission should come from identical zones within the nebula. The 4 He + and 4 He ++ 
emission can be measured using recombination lines. Unfortunately there is no direct way of 
measuring the amount of neutral helium within the H II region. In some cases the spectral 
types of the ionizing OB stars can be identified from optical photometry and spectroscopy. 
The helium ionization properties ar e then derived usi ng stellar atmosphere models to calcu- 
late the expected escaping flux (e.g.. lVacca et al.lll996l ). toget her with phot oionization models 



of the nebula to determine the He ionization structure (e.g., lRubinlll984h . At optical wave- 
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lengths spectral transitions of several other atomic species are also observed to help constrain 
these models. The most detailed analyzes have been made for the Orion nebula in order 



to determine the total 4 He/H abundance ratio (IMathis &: Rose 



Rubin et al. 1119911 IPogge et al. 111992 : Esteban et al 



1998; 



1991 



Baldwin et al. 1991 



Blagrave et al.l 120071 ) . Never- 



theless, the limiting factor in determining accurate total helium abundances in Orion is in 
converting 4 He + /H + to 4 He/H. 

Most H II regions are ionized by several stars and may thus have a complex geometry. 
Because of this little direct information may be known about the rad iation field. A va r iety o f 
diagnostics have been developed to probe the ionization structure. IVflchez &: Pagell (119881 ) 
developed a radiation softness parameter based on the fine structure lines of O and S ions 



that i s not v ery sensitive to chem ical composition (see also lShields fc Searle Ill978l ; IMathis 



19821 . Il985l ). I Armour et al.l (119991 ) use infrared transitions of Ne and Ar ions to determine 



the ionization structure. In general H II regions ionized by a hard radiation field produce 
the most accurate 4 He/H abundance ratios since all of the helium within the H II regions 
will be ionized. For example, metal poor blue compact galaxies with H II regions ionized 



of these objects has been extensively studied (teallantvne et al 



Gruenwald et al. 2002; Sauer fe Jedamzik 200 



Peimbert et al 



2000: 


Vieeas et a . 


2000 


2002: 


Izotov et al. 


2007) 



In some cases the total 4 He/H abundance ratio will be less than 4 He + /H + by a few percent, 
significant for cosmological implications, since the He II zone is larger than the H II zone. If 
clumping exits in these H II regions, how ever, the determined values will be underestimated 
by a few percent (IMathis fc Woodll2005l ). In the Galaxy the high excitation H II regions M17 
and S206 have been used instead of Orion to determine the total 4 He/ H abundance ratio 
(1Peimbertl[l~99~3l : lEsteban et al1ll999l : beharveng et alJbood : lBalseril2006h . 

Many of our sources are located throughout the Galactic disk and are totally obscured 
at optical wavelengths by dust. Furthermore, most of the optically visible H II regions in 
our sample are low emission measure H II regions in the outer Galaxy. The low EM renders 
useless the Orion-type detailed models since the spectral diagnostics used to constrain them 
are not available. Fortunately, unlike many of the light ele ments, 3 He/H abundan ces accurate 
to ~ 10% can yield important astrophysical conclusions (IWilson fc Roodlll994l ) . Because 
of this we adopt here a simple ionization structure model that is constrained primarily by 
H and 4 He radio recombination line observations. We then use numerical models to assess 
the accuracy of this approach for the H II regions in our sample. Our goal is to identify a 
subset of sources in our nebular sample wherein we can derive 3 He/H abundances accurate 
to ~ 10%. 
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1.3. 



3 He Ionization Correction 



We seek an ionization correction factor to derive the 3 He/H abundance from the 3 He + /H + 
abundance gotten from the 3 He + observations. We define the ionization correction factor, 
Ki, to be y 3 = Ki y£, where y is the He/H abundance ratio by number, the subscript denotes 
the isotope, and the superscript is the ionization state. 2 The ionization correction factor is 
determined by using the H and 4 He radio recombination line observations. We assume that 
the 3 He + emission traces the 4 He + emission, that the amount of neutral and doubly ionized 
helium is negligible, and a canonical value for the total 4 He/H abundance, y±. The 3 He/H 
abundance ratio is then given by 



Here £/4gal is the actual 4 He/H abundance by number in the Milky Way. There are, of course, 
some problems with this simple model. 3 He + and 4 He + are measured using transitions that 
are sensitive in different ways to density structure within the nebula. The 3 He + is observed 
using a collisionally excited hyperfine transition that is sensitive to the column density and 
thus proportional to the electron density, n e . The 4 He + is observed using recombination 
transitions that are proportional to the emission measure or n 2 e . Thus if the H II region has 
large density fluctuations the different transitions may probe significantly different material. 
For example, consider a simple two-component nebula with a small, very dense core and 
a larger, diffuse halo. Under certain physical conditions the halo will dominate the 3 He + 
emission while the core will dominate the 4 He + emission. If this is coupled with ionization 
structure between the core and halo then the simple formula in equation (1) will be incorrect. 

It has also proven difficult to measure accurately the 4 He/H abundance, y±, in the 
Galaxy. 4 He cannot be directly measured in the Sun and must be inferred from theoretical 
stellar evolution models and helioseismology. Measurements of 4 He in H II regions using 
recombination lines require an ionization correction as discussed above. A canonical value of 
2/4 = J/4GAL = 0.1 is adopted here although there is evidence that this often cited value may 
be too high (see §5}. For this reason we have parameterized equation (1) with the 2/4GAL 
factor so that our ionization corrections can be easily scaled should an accurate Milky Way 
?/4 value be derived in the future. 



2 Thus the notation used throughout is: y% = 3 He/H; y% = 3 He + /H + ; y± = 4 He/H; and = 4 He + /H+. 




(1) 
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1.4. Radio Recombination Line Observations 

Here we use high signal-to-noise ratio recombination line spectra taken at two different 
spatial resolutions to probe for any significant ionization structure. For our H II region sample 
we have obtained H and 4 He RRL data simultaneously for several different transitions near 
8 GHz with a spatial resolution of 3'.5 (see Paper I). These observations are briefly reviewed 
here in § 12.11 Because they sample the same H II region zone as the 3 He + spectra, these 
data are primarily used to determine y\ in equation (1). Most 4 He RRL data are from 
single-dish telescopes which typically have spatial resolutions > 1'. Measuring accurate 
y^ abundance ratios with single-dish telescopes is difficult because non-random frequency 
structure is produced in the instrumental baselines ow ing to reflections from various parts of 



the telescope structure (see lLockman fc Brown! Il982l ). Improvements in receiver technology 



have provided better stability and enhanced signal-to-noise. This, together with a better 
understanding of the instrumental effects, has enabled more accurate measurements of y\ 
(Peimbert et al. 1992a; Paper I). 

Using these improved techniques we made additional observations at 18 GHz with a 
spatial resolution of 1'5 for a subset of our H II region sample (§ 12.21) . The smaller beam of 
the 18 GHz RRLs typically probes the more compact, dense components, whereas the larger 
3f5 beam of the 8 GHz RRLs is more sensitive to the extended, diffuse material. Radio 
interferometers can measure the 4 He + distribution within H II regions at spatial resolutions 
< V . Although interferometers do not detect the diffuse emission because of the zero- 
spacing flux problem, this information can be provided by single-dish measurements. To 
our knowledge only five H II regions in our sample have been observed in H and 4 He RRL 
emission with radio interferometers: W3, SgrB2, W43, W49, and W51. 

Since the entire nebula can be imaged, the best optical data for this study are obser- 



vations made with Fabry-Perot spectrophotometers. ICaplan et al.l (120001 ) derived oxygen 
and helium abundances for 34 H II regions using the ESOP Fabry-Perot instrument. For 
many objects the entire nebula was probed ( the largest diaph r agm w as 4' 27"). The 4 He 



abundances from this survey are discussed by iDeharveng et al.l (120001 ). There are four H II 



regions that are in our sample: S206, S209, S212, and S252. 

Although the 8 GHz RRLs are primarily used to determine y\ the other radio and optical 
4 He + data can not only provide an independent confirmation of the 4 He + /H + abundance, 
but also indicate whether the assumptions in equation (1) are valid. In §[3] numerical models 
are used to explore the physical conditions under which the assumptions in equation (1) 
begin to fail. This information is then used in §H] to determine a value for the ionization 
correction factor, Ki, for each source. A discussion of the results is in §[5] and a summary of 
the paper in §0 
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2. 4 He+ OBSERVATIONS 

2.1. X-band (8 GHz) Radio Recombination Lines 

The X-band RRL observations are discussed in §2 of Paper I. The H91a and He91o 
transitions were observed in every H II region. Additional, higher order, transitions were 
also observed at nearby frequencies with essentially the same spatial resolution (3f5). These 
data were used to constrain the nebular models in Paper II. Here we only use the higher 
quality 91a and 114/? transitions of H and 4 He. 

Calibration details are discussed in §2.1 of Paper I. We used the planetary nebula 
NGC 7027 to calibrate the flux density scale for each observing epoch. The absolute calibra- 
tion is judged to be within ~ 10%; the internal consistency is ~ 5%. From these observations 
we use only the ratios of the line areas to calculate 4 He + /H + ; the calibration is therefore not 
a significant source of uncertainty. 

A much larger source of uncertainty are the standing waves produced by reflections 
from various parts of the telescope structure (see §2.2 in Paper I and references therein). To 
minimize the amplitude of these standing waves, observations are made with focus offsets 
alternating between ±Ao/8 every 60s. Cancellation of the standing waves is, however, not 
perfect and the spectral baselines must be modeled. A 12 th order polynomial baseline is used 
to model this non-random frequency structure. The necessity for usi ng a high-order bas eline 



and the selection of the appropriate order to be used is discussed in iBalser et al.l (119941 ) and 



Paper I. The results of Gaussian fits to each transition are given in Table 4 of Paper I. 



2.2. K-band (18 GHz) Radio Recombination Lines 

Observations of H and 4 He RRLs at K-band were made with the NRAO 140 Foot 
telescope during the periods: 1995 November, 1995 December, and 1996 January. The 140 
Foot telescope has a half-power beam-width of 1'.5 at 18 GHz. The radiometer consisted of 
a dual circularly polarized HEMT receiver located at the Cassegrain focus of the 140 Foot 
telescope. The system temperature on cold sky under good weather conditions was ~ 50 K. 
The spectrometer consisted of the NRAO model IV 1024 channel autocorrelator (AC). The 
AC sampled four independent 20 MHz wide frequency bands ( "quadrants" ) of 256 channels 
each. The velocity resolution at 18 GHz is 1.3 km sec -1 . The 20 MHz bandwidth allowed both 
the H and 4 He transitions to be observed simultaneously within each quadrant. For most of 
the observing the AC was configured with the H70a/He70a and the H88/?/He88/? RRLs in 
two quadrants each, simultaneously sampling orthogonal polarizations in order to maximize 
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our spectral sensitivity. In some of the brighter objects the IOO7 and 1105 transitions of H 
and 4 He were observed. 

Spectra were taken using the total power observing mode. First a reference spectrum 
(OFF) was taken at a position offset 6 minutes in right ascension from the nebula and then 
the nebula itself (ON) was observed. Each position had an integration time of 6 minutes. 
Local pointing corrections were determined every hour. After each pointing a continuum 
cross scan, in right ascension and declination, was made on the H II region to measure the 
radio continuum emission. 

The radiometer system temperature was measured using noise tubes calibrated by con- 
necting the receiver to matched hot and cold loads. An ambient temperature absorber was 
used for the hot load and the zenith sky for the cold load. Since this procedure has many 
uncertainties and is not always performed before each observing session we have re-calibrated 
the system temperature by using the observed continuum intensity of the unresolved plane- 
tary nebula NGC 7027. NGC 7027 was observed several times during eac h observing sessio n 



during good weather near transit; it has a flux density of 6 Jy at 18 GHz (jPeng et al. Il2000l ). 



There are two effects that can modify the calibration as a function of time. First, at 
these frequencies the aperture efficiency of the 140 Foot telescope changes with elevation as 
gravity distorts the primary reflector; an elevation gain correction must therefore be made. 
We measured the H II region continuum intensity every hour for each source in order to 
determine the appropriate gain correction. Second, the atmospheric opacity is significant 
and depends on the water vapor content. We made a tipping scan several times during each 
observing period to measure the opacity. Since we are interested in calculating the 4 He + /H + 
abundance ratio, which depends on the ratio of the line intensities, the calibration effectively 
cancels. All intensities quoted here are therefore antenna temperatures calibrated using only 
the NGC 7027 data. 

We used the techniques discussed in §2.2 of Paper I to model the K-band spectral 
baselines, except that a 6 th order polynomial was used instead of a 12 th order. Although 
the AC configuration was the same, yielding a 20 MHz bandwidth with 256 channels, only 
~ 12 MHz is required to cover the frequency range between the H and 4 He lines. We have 
therefore reduced the total number of terms in the polynomial baseline model. 



2.3. The 4 He+ / H+ Abundance Ratio 

Example RRL spectra are shown in Figured! Only 10 objects from the original sample of 
21 were observed at K-band. The data were smoothed to a velocity resolution of 8.1 km sec -1 . 
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The vertical lines flag the expected locations of the H, He, and C lines for the specified 
RRL transition (91a or 70a). Table [1] summarizes the results of Gaussian fits to the K- 
band spectra. Listed are the source name, the RRL transition, the peak intensity and full- 
width half-maximum (FWHM) line-width together with their associated errors, the total 
integration time, the resulting r.m.s. noise of the spectral baseline, and the quality factor 
of the line. The quality factor is determined using the same criteria as in Paper I. These 
include the signal-to-noise ratio, the structure of the baseline, and the crowding of spectral 
lines. We do not include any calibration error in assessing the quality factor since we only 
use line ratios in this paper. 

The 4 He + /H + , abundance ratios calculated from the RRL data using the line areas 
are summarized in Figure [2] which plots the y\ abundance ratio as a function of the principal 
quantum number, n, of the RRL transition. The 70a, 88/3, 91a, and 114/3 transitions are 
shown. The Gaussian line parameters in Table 4 of Paper I and in Table 1 here are used for 
the X-band and K-band transitions, respectively. The errors are determined by propagating 
the uncertainties in the Gaussian fit to the line intensity and width. 

For some nebulae there are real differences among the y\ ratios derived from the various 
RRL transitions. This is to be expected, of course, since the K-band spectra probe, on 
average, the denser, more compact nebular interior whereas the X-band spectra measure 
a much larger region, perhaps including a low density halo. These statistically significant 
differences are clearly seen in Figure [3] which compares in various ways the y\ ratios derived 
from the H91a (diamonds) and H70a (circles) spectra for the 10 K-band sources (see §EJ). 

3. SYNTHETIC NEBULAR MODELS 

We use the formalism defined in §1.31 together with numerical models for the nebulae to 
explore the range of validity of this approach to the ionization correction. The H II regions in 
our sample span a wide range of physical conditions. At one extreme there are large, complex 
H II regions such as W49 that are ionized by many OB-type stars and have multiple compact 
components in their density structure. At the other extreme are H II regions such as S206, 
which are smaller, much less complex sources, that are probably ionized by a single OB-type 
star. The simple ionization structure correction model discussed in § 11.31 could well not be 
very accurate for sources like W49. It is important, however, to understand under what 
conditions this model breaks down and what impact this has on deriving 3 He/H abundance 
ratios. 

Paper II explored the effect of source density structure on the derivation of the 3 He/H 
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abundance. We used RRL and continuum data to constrain the models. Both RRL and con- 
tinuum emission intensities are proportional to n\ and so they trace the densest components 
of the nebulae. We found that our sources could be grouped into two broad categories. In one 
class were "complex," W49-like sources that required models with complex density structure 
and non-LTE radiative transfer in order to reproduce the observations. In the other class 
were "simple," S206-like sources whose observations could be reproduced by models with 
LTE radiative transfer and homogeneous density structure. 

We adopt a simple core-halo model for the nebular density and ionization structure. 
From the models we calculate simulated observations and assess the effect of the model 
structure on our ability to recover accurately the 3 He/H abundance using our simple ioniza- 
tion structure correction. We use a computer program called NEBULA to model the H II 
region (see Balser 1995 and Paper II for details about NEBULA). The model nebula can have 
arbitrary density, temperature, and ionization structure. The gas consists of only hydrogen 
and helium. We calculate the radiative transfer through the model nebula for the radio con- 
tinuum, recombination line, and 3 He + hyperfme line emission. Since we shall be comparing 
these synthetic spectral lines to the §|2]data taken at 8 and 18 GHz, we calculate the model 
antenna temperature for each species by convolving the model sky brightness temperature 
distribution with the telescope beam and multiplying by the telescope main beam efficiency 
The telescope beam is assumed to be a Gaussian with a half-power beam-width (HPBW) of 
1^5 or 3'5. These are the HPBW sizes of the 140 Foot telescope beam at 18 GHz and 8 GHz, 
respectively. 

All model nebulae are located at a distance of 8 kpc which is roughly the average dis- 
tance to our 3 He sources. The models are comprised of nested, spherical, homogeneous, 
isothermal components. For simplicity, all model components have a fixed, 8,000 K , electron 
temperature that is typical for Galactic H II regions. For all models it is assumed that the 
total abundance ratios are y 3 = 1 x 10~ 5 and y± = 0.10. Table [2] summarizes the model 
parameters used as input for the NEBULA program. 

Model A is a one-component, low density, n e = 100 cm -3 , nebula with an angular size of 
3'.5 (8.14 pc). All other models are two-component (core-halo) nebulae which consist of a l'.b 
(3.49 pc) core embedded in a larger 3'.5 nebula. Listed for each component is the angular size, 
the electron density, the fraction of the component H II mass relative to the total mass, the 
Hz abundance ratio, and the y~X abundance ratio. Altogether, the models span a core-halo 
density contrast ratio from 1:1 to 5:1, and a core mass fraction from 8% to 30%. That is, all 
models have halos that contain more mass than the cores by a factor of either ~ 2 or ~ 10. 
These parameters span the range of the Galactic H II regi on physical properties derived from 
single-dish RRL observations (e.g., Quireza et al. 2006al Jbh . 
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The NEBULA code calculates the following observable quantities: the continuum an- 
tenna temperature at 8.665 GHz and the antenna temperatures of the 3 He + , H91a, He91a, 
H70a, and He70a transitions. Properties of the synthetic observations that result from the 
Table [2] models are summarized in Table [31 Listed are the continuum and line antenna tem- 
perature for the 3 He + transition, the y£ ( 3 He + /H + ) abundance ratio, the y\ ( 4 He + /H + ) 
abundance ratio and the values of Ki and y$ calculated using the equation (pQ) Ansatz. 
All of these quantities are determined for X-band (8 GHz). The last column lists the y± 
model values at K-band (18 GHz). The y\ values are calculated from the H and 4 He line 
parameters produced by NEBULA; they are the ratio of the line areas. The yt, value is 
derived from the NEBULA models assuming a single, homogeneous, fully ionized sphere 
where yt, oc Tl Av/T^ 2 . Th is analytical expression approximates the numerical calculation 



to within a few percent (see iBalserl Il995l ) . This is our standard way of deriving the 3 He/H 
abundance ratio (Paper II). Here we apply this analysis to the Table [2] nebular models and 
use the formalism developed in §1.31 to derive the ionization correction factor, K{. 

Model Al is a uniform nebula with no density or ionization structure. The NEBULA 
model abundances in Table El yt and yt, exactly match the input abundances, y% = 1 x 1CT 5 
and ?/4 = 0.10. This of course must be the case since Model Al is constructed to match our 
analysis assumptions. It is a single, homogeneous, fully ionized sphere and the NEBULA 
Model Al gives an ionization correction of hZi — 1 clS it should. 

In Model A2 we simulate a one-component nebula that has 30% of all the helium neutral. 
For this model the NEBULA 4 He RRL synthetic spectra together with equation (CQ) give an 
ionization correction of = 1.4 which in turn produces y$ = 1.0 x 10 -5 , the model input 
abundance. Thus, for a one-component model nebula the observables are sufficient to derive 
ionization corrections that can accurately recover the input structure. 

Model B has a core-halo density structure where the core component has a higher density 
by a factor of 5, but the halo has a factor of ~ 2.3 times more mass. Both core and halo 
are fully ionized: Model B has no ionization structure so we cannot derive any correction 
factor. The NEBULA Model B yt abundance thus underestimates the true yt abundance 
by a factor of 1.3. In Paper II we discuss density structure corrections for the 3 He + /H + 
abundance ratio derivation. There we model nebular density structure and estimate the 
density correction factor for each H II region in our sample. Applying these techniques to 
Model B recovers, as expected, the input yt abundance. 

Models C and D are constant density nebulae with core-halo ionization structure. The 
halo mass in both models is ~ 11.5 times that of the core. Model C has a fully ionized core 
and a halo that is under-ionized by 50%. Model D reverses this ionization structure. That 
these models have ionization structure is revealed by the different NEBULA values for the X 
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and K-band abundance ratios (see Table [3]). Since most of the mass resides in these halos 
the NEBULA 3 He + /H + abundance ratio is much less for Model C with its under-ionized halo 
(yt = 0.57) than for model D (y^ = 0.96). The NEBULA results give ionization corrections 
of Ki= 1.82 and 1.05 for models C and D, respectively. In both cases Equation [T] yields y 3 
= 1.0 x 10~ 5 , the model input value. For these two models our simple ionization correction 
reproduces the correct 3 He/H abundance ratio. 

Models E and F combine both density and ionization structure in a core/halo geometry. 
Both have high density cores (5:1 density contrast) and halo masses that are ~ 2.3 times 
that of the cores. Model E has a fully ionized core and a halo that is under-ionized by 50%. 
Model F reverses this ionization structure. Once again the different NEBULA values for the 
X and K-band y\ abundance ratios reveal the presence of nebular ionization structure. For 
model E the NEBULA y^ value is 0.53, roughly a factor of two too low. Yet the NEBULA 
y^ abundance ratio implies only a small ionization correction factor, = 1.18, which gives 
a corrected y% abundance of only 6.2 x 10~ 6 , 62% of the input value. The situation is much 
better for Model F with its fully ionized halo. The NEBULA values for y± and y% together 
give an ionization correction factor of Ki = 1.56 and a y% abundance that is within 2% of the 
input value. 

These models show that the simple ionization correction given in Equation [1] will be 
reasonably accurate for nebulae with core/halo geometries. Our correction begins to break 
down whenever the low density component contains most of the mass and is under-ionized. 
This is because the 3 He + and RRL emission arises from gas located in different zones of the 
nebula and these regions do not have the same He ionization structure. 

In sum, the simple ionization correction given by equation ([TJ) provides a reasonable value 
for the total y^ abundance ratio except for the Model E H II region that has an under-ionized 
halo containing a significant fraction of the total mass. Clearly more complicated density and 
ionization structures exist in real H II regions. The simple models discussed here nevertheless 
illustrate the main issues. The simple ionization correction breaks down when ionization 
structure exists between significantly different density components. True variations of y^ 
and 1/3 structure within the H II region will complicate the analysis. Nonetheless, because 
of our large X-band beam that in many cases covers the entire nebula, emission from dense, 
compact components tends to be diluted. We thus derive an average value of y% for each 
object. 
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4. IONIZATION CORRECTION FACTOR 

Here we assess the ionization structure correction for our 3 He sample of 21 Galactic 
H II regions. In Paper II we classify nebulae into two categories: (i) "simple" H II regions 
dominated by low-density, extended structure that is well fit by a single-component uniform 
density sphere model; and (ii) "complex" H II regions that have multiple high density com- 
ponents and thus require a more sophisticated model. The simple H II regions are: G133.8, 
S206, S209, S212, S228, S235, S252, S311, S76, S90, S156, NGC 7538, and S162; the complex 
H II regions are: W3, SgrB2, Gl.l, M17S, M17N, W43, W49, and W51. Simple H II regions 
characteristically have low emission measure. For complex H II regions that have significant 
ionization structure we will underestimate Ki as discussed in §31 We therefore can only derive 
a lower limit to the 3 He/H abundance for these nebulae. 

We used various techniques to assess the excitation state of the nebulae. Table 0] sum- 
marizes the excitation properties for our H II region sample. Listed for the radio data are the 
4 He + /H + abundance ratio, y^, derived using only the H91a and He91a RRLs, the H-ionizing 
luminosity, Al, determined from the radio continuum data in Paper I, and the spectral type 
of a single ZAMS star that can produce t he observed luminos ity. Spectral types a re derived 



from stellar atmosphere models: both the IVacca et al.l (119961 ) and iPanagial (119731 . m paren- 
theses) models are listed. Similar results from optical data are taken from the literature. The 
spectral types are based on photometry and spectroscopy; the earliest type star thought to 
ionize each nebula is listed in Table HI Overall, the radio and optical spectral classifications 
are in good agreement. 

We plot the ionization state, y£, of the 14 nebulae for which we have sufficient infor- 
mation as a function of the log of the H-ionizing luminosity, Al, in Figure HI There is no 
discernible trend connecting y\ and Al. Except for their systematically lower Al values 
(recall they are low emission measure nebulae), the simple sources span the same y\ range 
as do the complex nebulae. Values of Al determined from the radio continuum include the 
contribution of H-ionizing photons from all of the early-type stars within the H II region 
complex. Therefore low emission measure does not necessarily imply that the radiation field 
is soft. For example the radiation field of the low emission measure Rosette nebula is domi- 
nated by 04 stars compared to the rather weak ionization field of the high emission measure 
Orion nebula. 

Below we discuss each H II region separately in more detail. An ionization correction 
factor, Ki, is determined for each object. For consistency we only use data as of 1996 March 
(cf., Papers I and II). The k% values can be multiplied b y the 3 He + /H + ab undance ratios 



in Table 5 of Paper II to yield the 3 He/H abundance. iRood et al.l (120071 ) derive 3 He/H 



abundances for all the 3 He nebulae using the entire 1982-1999 140 Foot dataset. 
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W3 — Although W3 is optically obscured it has been extensively studied at radio 
wavelengths. Our higher spatial resolution K-band data give y\ = 0.081 whereas the X-band 
data yield y\ = 0.076. Both data sets are of very high quality. One interpretation of these 
results is that the more compact components of W3 — in particular W3A — have a higher 
excitation than the outer regions probed by the larger X-band beam. This is sometimes 
called the geometric effect where yX ratios increase with higher resolution observations. 

Interferometers can probe the RRL emission of the nebular gas at high spatial resolution. 
Extant observations include: 76a and 110a data fr om the Very Large Array (VLA) and the 
Westerbor k Synthesis Radio Te l escope (WSRT) (Roelfsema fc Gossl 1991 ) and 92a from 
the VLA (IRoelfsema et al.lll992l ; lAdler et al.lll996l ). The 4 He + /H + abundance ratio varies 
considerably among the different components (0.06 < y\ < 0.20). Values of y\ higher than 
0.10 are thought to be due to local enrichment from a nearby evolved object. Alternatively. 



radia tive transfer effects may cause one to overestimate the true y^ abundance (iGulyaev et al. 



19971 ) . High yt values in other H II regions are associated with evolved early-type objects 



[e.g.. iBalser et al.ll200ll ). W3 is rather unusual in our sample since the total mass of ionized 
gas is small (20 M Q ; paper I) and therefore mass loss from local objects can significantly 
pollute the surrounding material. We nevertheless adopt y\ = 0.076 as a reasonable value for 
most of the mass in W3. Because W3 has complex density structure and there is uncertainty 
about the ionization structure, we adopt an ionization correction of Ki > 1.32. 

G133.8 — Also known as W3N, G133.8 is the northern most region of the W3/W4 com- 
plex and has a slightly more negative velocity (by 7 km sec -1 ) than the other components of 
W3. Infrared and radio continuum observations of G133.8 indicate that this region is ionized 



by a l uminosity equivalent t o an 07 star (Table HI ISchraml fc Mezgerlll969l ; iThronson et al. 
19841 ; ICarpenter et al.l 120001 ). The X-band RRLs measure y\ = 0.082, roughly consistent 



with the ionization properties predicted from the continuum observations. We adopt an 
ionization correction of k^ = 1.22. 

S206 — This outer Galaxy H II region is thought to be ionized by a single star with 
spectral type 04-06. Fabry-Perot spec trophotometer data sh ow that S206 contains essen- 
tially no neutral helium with yf = 0.10 (IDeharveng et al.ll2000i ). Our X-band measurements 
give a slightly lower value of y\ — 0.092. Rece nt high sensitivity RRL data with the Green 
Bank Telescope yield y\ = 0.085 (jBalserl 120061 ) . Since we are correcting 140 Foot 3 He data, 
however, we adopt Ki = 1.09. 

S209 — This H II region has the largest Galactocentric distance, i? ga i = 16.9 kpc, in our 
entire sample and t hus is an important o bject. Spectrophotometry reveals an 09 and two Bl 
stars exciting S209 fjChini &: Wink! 1 19841 ); this suggests that S209 m ight contain some neutra l 
helium. Optical 4 He observations show that y\ ranges from 0.084 (jVflchez fc Estebanlll996l ) 
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to 0.12 (jDeharveng et al.ll2000i ). The He91a transition gives y\ 
0.083 producing an ionization correction of Ki = 1.20. 



0.083. We adopt y% 



S212 — 3 He + was not detected toward S212 and the 4 He data are of poor quality. S212 
is nonetheless a potentially important object with R ga \ = 14.2 kpc. The nebula is ionized 
by an early-type star, 05.5-08. Optical observations measure y\ = 0.088-0.10; most of 
the helium should be ionized. We adopt y\ = 0.088 which gives an ionization correction of 
Ki = 1.14. 

S228 — 3 He + was not detected toward S228 and the He91a line is only a probable 
detection. As is the case for S212, this outer Galaxy nebula is potentially a good 3 He 
candidate. Based on both radio and optical data, however, this region is ionized by a late- 
type O-star; the helium is probably not fully ionized within the H II region. Because 3 He + 
and 4 He + have not been detected we have a direct measure of neither the 3 He/H abundance 
nor the ionization correction. 

S235, S76, S90 — These physically small, shell-like nebulae were observed because they 
have a morphology that is similar to that of W3 A. The anomalously high 4 He/H abundance 
ratio in W3A might be caused by local enrichment from nearby, evol ved objects. This might 
also explain the high W3 3 He/H abundance ratio (lOlive et al.lll995l ). 



3 He + was not detected in either S235 or S76. Although we do not have good quality 
4 He + data it is clear that the helium in these nebulae is significantly under-ionized (Paper 
I). Both radio and optical data show that these two H II regions are ionized by late O-type 
or early B-type stars. Given these soft radiation fields even a high 3 He/H abundance ratio 
would be difficult to detect. 

In S9 0, both 3 He + and 4 He 4 " are detected. Although an 09.5 star is optically identified 



near S90, ICrampton et al.l (119781 ) argue that this star does not excite the nebula. The 
radio continuum data suggest an earlier type object, 08-06.5, roughly consistent with the 
observed value of y\ = 0.070. We adopt an ionization correction factor of Ki = 1.43. 

S252 — Also known as NGC 2175 or W13, this diffuse H II region extends over 25' 
on the sky and contains several compact c omponents. T he nebula is ionized mainly by an 
06-06.5 star centered on the H II region (lHaikalalll995l). Our observ ed position lies on a 
bright rim located at the western edge of the nebula (IFelli et al.lll977l ). Because of the low 



emission measures in S252 there is only a probable dete c tion o f He91a (Paper I). 4 He + is, 
however, detected at optical wavelengths. IShaver et al.l (119831 ) measured y^ = 0.051 near 



our p osition. Towards the center of the H H re gion this value increases with y^ = 0.091- 
0.10 (IShaver et al.lll983l ; iDeharveng et al.l |2000| ) . We adopt an ionization correction factor 
of Ki = 1.96. 
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S311 — S 311 is thought to be ionized by an 06-06.5 star located north of the radio 



emission peak ( Persi et al 



type geometry ([Albert et al 



987) 



The asymmetric density distribution suggests a blister 
19861 ). He91a emission is detected: yX = 0.071. We adopt an 
ionization correction factor of Ki = 1.41. 

SgrB2 — SgrB2 is the most complex source in our sample. It is important because it 
is located near the Galactic Center (G C). The ionized gas is observed at high spatial res- 
olution with radio interferometers fe.g . . iMartin fc Downed Il972l ; iBenson fc Johnston! Il984l ; 

1995h. Ther e are at least 57 separate, compact com- 



Gaume fc Claussen 1990; Gaume et al. 



ponents associated with SgrB2 (IDePree et al.lll996l ). 



The 4 He + /H + abundance in SgrB2 is known to be anomalous (see lLockman & Brown 



19821 . and references within). The radio continuum emission implies a very high excitation 
but the 4 He + /H + ratio is at times reported to be below 0.06 (e.g., TableH]). In the literature 
the SgrB2 4 He/H abundance is reported as underabundant, overabundant, or normal, de- 
pending on various effects such as the geometric effect, selective photo- absorption by dust, or 
the external maser effect. These uncertainties are compounded by complex spectral baselines 
produced by reflections from the telescope structure for the single-dish observations. 



Some insight can, however, be gleaned from high resolution inter: 
of RRL emissio n. Sgr B2 was observed with the VLA in 76a emission 



erometric observations 



Roc 



sema et al.ll 19871) ; 

110a emission (iMehringer et al.lll993l ); and 66a emission (IDePree et al.lll995l . Il990 ). The 
4 He + /H + rati o varies from 0.05 to 0.17 among the yariou s compact components; the average 



value is 0.10 (IRoelfsema et al.lll987l ; IDePree et al.lll996l ). The component F bright cluster 
of so urces is the sole extremely low 4 He + /H + ratio, y\ < 0.05, component (IDePree et al. 
1996h . 



Most of the radio continuum flux d ensity at X-band arises from the compact compo- 
nents towards SgrB2 (IBalser et al.lll995l ). The SgrB2 F components have the highest peak 
emission measures. Since the RRL emission is proportional to emission measure, the 4 He + 
RRL emission is strongly influenced by the F components. It is therefore not surprising that, 
depending on frequency and spatial resolution, different yX values are measured in SgrB2. 
Because toward Sgr B2 there is ionization structure coupled with very complex density struc- 
ture within the 3'5 beam of the 140 Foot telescope, using equation (1) is inappropriate. For 
this source we can only adopt a limit to the 3 He/H abundance ratio. 

Gl.l — Better known as SgrD, this H II region is located only 1° from the GC but is 
much less complex than SgrB2. Although we have placed its location at the G C there is 



evidence suggesting that it might be located outside the molecular nuclear disk (ILisI Il991 



Mehringer et al.lll998t iBlum &; Damineh1ll999l ). High resolution radio continuum images re- 
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veali several compact components within a d iffuse . exten ded halo (jLisztlll992l ; lMehringer et al. 
19981 ). A single 07 star can ionize the core (ILisztlll992l). This is con sistent with narrow band 
imaging at 2.17/im, Br7, and 2.06/xm, and He I (IBlum et al.lll999l ). If the diffuse gas is in- 
cluded in the analysis an 04 star is required to account for the excitation (Table Hj). The 
4 He + /H + ratio of 0.064 implies the presence o f some neutral heli um in the nebula. This is 
also consistent with the infrared observations of IBlum et all (Il999l ). Initially Gl.l was classi- 
fied as complex because of scant data; additional RRL data supports a simple homogeneous 
model. We therefore adopt the ionization correction factor of = 1.56. 

Ml 7 — The Omega nebula, also known as M17, S45, NGC 6618, and W28, is a well 
studied inner Galaxy H II region. The nebula spans over 10' in angular size and consists 
of two main components: M17 north (M17N) and M17 south (M17S). The Omega nebula 
i s complex: high resolution radio continuum images reveal this density structure directly 
( jFelli et al.lll984j ). We were unable to model the density structure in Paper II. The complex 
velocity field produces non-Gaussian profiles in both M17N and M17S (Paper I; also see 
Joncas & Roy 1986). 

It is nonetheless clear that most of the helium is ionized. High extinction makes optical 
identification of t he ionizing stars difficu lt, but infrared spectrophotometry shows at least 
five 03-06 stars (IHanson fc Contil Il995l ) . The large number of early O-type stars and the 
high 4 He + /H + abundance ratios observed at both radio and optical wavelengths suggest that 
Ki ~ 1. Using the He91a transitions we adopt an ionization correction factor of = 1.09 
and 1.06 for M17S and M17N, respectively. 

W43 — W43 is a large, optically obscured H II region locat ed only 4.6 kpc from t he GC. 
Because H II regions are sparse between the GC and 4.5 kpc (ILockman et al.lll996l ). W43 
probes the 3 He/H abundance near the outer edge of this H H region desert. W43 contains 



low density, diffus e gas and many compact components (e.g., ISubrahmanyan fc Gosslll996 



Balser et all l200ll ). Near infrared observations show a star cluster at the c enter of the 



H II region consisting of a WN7 star and two O-type stars (IBlum et al.lll999l ). The total 
H-ionizing luminosity is considerable, requiring at least an 04 star for the entire region 
(Table i}. 

We studied the helium ionization properties of both the diffu se and compact gas in 
W43 using the 140 Foot telescope and the VLA (IBalser et al.ll200ll ). Using RRL emission 
we imaged the 4 He + /H + abundance ratio for both the compact and diffuse components. The 
average 4 He + /H + ratio of the entire dataset is (y^) = 0.077 ± 0.01. We find no significant 
variations of the 4 He + /H + abundance, even for positions observed at the edge of the nebula 
with the 140 Foot. Thus although W43 is a complex nebula, the ionization properties do not 
appear to vary within this object so we can derive an accurate ionization correction factor. 
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We adopt = 1.30 based on the average y\ abundance. 



W49 — W49 is one of the most luminous H II regions in the Galaxy, consisting of a 
ring of ult ra compact H II r egions (UC H II regions), each containing at least one O-type 



star (e.g., IWelch et al.l 119871 ). T hese UC H II region s make up only one-third of the total 



H-ionizing luminosity observed (jConti fc Bfuml 120021 ) . The He91a RRL gives a 4 He + /H + 
abundance ratio of 0.079 ± 0.0052, consistent with the K-band data. Higher resolution VLA 
H and He RRL images find y t ranging from 0.05 ± 0.02 to 0.18 ± 0.06 in 13 objects with an 
average of (yf) = 0.11 ±0.01 (IDePree et al.lll997l ). This significant ionization structure com- 
bined with considerable density structure makes any equation ([1]) type ionization correction 
inappropriate. We therefore adopt Ki> 1. 

W51 - - W51 is a bright, complex H II region with an angular extent of 2° on the 
sky. Radio continuum observations reveal several 04-06 stars associated with the main 



components (IMehringerl Il994j) . Mid- infra red imaging and spectroscopy suggests stars with 
spectral types ~ 09 (jOkamoto et al.ll200ll ). The X-band 4 He + RRL data give a yf abundance 
of 0.079 ± 0.0024. The He70a transition appears to be slightly lower with y\ ~ 0.07. 
High resolution VLA observations of 92a RRLs yield yt values ranging from 0.08 ± 0.02 to 
0.14 ± 0.06 with an average of (yt) = 0.096 ± 0.015 J]VIehringerlll994h . Although most of 
the helium appears to be ionized, given the complexity of this object we adopt an ionization 
correction factor of Ki > 1. 

S156 — S156 is a compact H II region excited by an 06.5-type star. We have only an 
upper limit for the 3 He + /H + abundance. Although we have detec ted 4 He + at X-band th e 
accuracy is not very good. Optical observations yield yt = 0.081 (IDeharveng et al.ll2000l ). 
We adopt an ionization correction factor of Ki = 1.23. 

NGC 7538 — Al so known as S158, NGC 7538 is centrally diffuse with a bright rim to 
1977h 



the west (llsrae 



An 07 star l ocated at the center of the neb ula is responsible for 



the e xcitation (IDeharveng et al.lll979l ). Optical data yield yt = 0.10 (ILynds fc O'Neil. Jr. 



19861 ). a value higher than our RRL results. Both the X-band and K-band data give yt 
~ 0.08, although the Hell4/3 transition produces a slightly higher value of 0.10. We adopt 
an ionization correction factor based on the X-band data of k, = 1.19. 



Sit 



Sometimes called the Bubble nebula, S162 consists of a circularl y shaped shell 



with bright emiss i on to the north observed in Ha and the radio continuum (IBarlow et al. 



19761 ; llsrael et al.l Il973l . Paper I). The exciting star has been classified as an 06.5-type, 
consistent with the radio continuum estimates. We do not detect 3 He + emission. Although 
we do detect the He91a RRL, the quality factor is poor. We cannot make an ionization 
correction so therefore adopt Ki > 1. 
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DISCUSSION 



Our nebular ionization correction, Ki , is the simplest possible approach to the problem. 
It makes several key assumptions: (1) the 4 He abundance is 4 He/H = 0.10 for all Galactic 
H II regions; (2) all the He is singly ionized — there is neither neutral nor doubly ionized 
He within the nebulae; and (3) the 3 He + and 4 He + ions occupy identical volumes within 
the nebulae. Measuring the total helium abundance is difficult since there is no direct way 
to observe neutral helium. We do not expect any significant doubly ionized helium since 
the radiation field from Galactic O-type stars is not hard enough. For example, M17 is 
ionized by several stars classified between 03-06 (IHanson fc Contil Il995l ) and is one of the 
Galactic H II regions with the highest degree of ionizati on. Yet only a very low upper limit, 
3 He ++ /H+ < 8 x 10~ 5 , is measured for 4 He ++ in M17 JPeimbert et al.lll992bh . In fact, to 
our knowledge 4 He ++ has never been detected in any Galactic H II region. 

A canonical value for the Galactic 4 He abundance of t/ 4 = 0.1 is typically used in the 
literature. Most of the 4 He was produced during the era of primordial nucleosynthesis. 
Therefore the variation of 4 He in the Galaxy, caused by stellar and Galactic evolution, is 
expected to be small bec ause AY/AZ ~ 1, wh ere Y and Z are the helium and metal 
abundances by mass (e.g., IChiappini et al. Ibooj ). Direct measurement of the total 4 He/H 
abundance ratio in the Galaxy is difficult. The 4 He abundance cannot be determined using 
either the Solar photosphere or in meteorites. Measurements of 4 He in Galactic H II regions 
have to correct for ionization structure just as we are trying to do here for 3 He. 

Solar 4 He abundances are determ i ned from t h eoret ical stellar evolution models with y 4 
~ 0.1 (e.g.. lAnders &: Grevesse Ill989l ). iPeimbertl (119931 ) suggested that the high excitation 
H II region M17 was the best object to measure 4 He directly and they also determined 
U4 ~ 0.1. There is mounting evidence, however, that y 4 is less than 0.1 in the Galaxy. More 
recent calibration of the 4 He Solar abundances yield y\ values lowe r by as much as 10% 
( jGrevesse et al.lll996l . Il998l ; iBasu fc Anita 112004 ; lAsplund et al.ll2005l ). Recent estimates of 
y 4 in M17 that inclu de temperature fluctuation s reduce ?/ 4 by about 5%, consistent with 
improved RRL data (IPeimbert fc Peimbertl 120021 ; lOuireza et al.ll2006al ). Furthermore, RRL 
observations of the lower metallicity Galactic H II region S206 , which is also expected to 
contain little or no neutral helium, give |/ 4 = 0.085 (IBalserl 120061 ) . 



Nevertheless, we adopt y 4 = ?/ 4 gal = 0.1 in equation (1) for our ionization correction 
factor analysis. We judge that for our simple sources the ionization correction should be good 
to about 10%. At this level of accuracy all of the additional ionization effects considered by 
4 He abundance analyses (e.g.. Ilzotov et al.ll2007l ) can be neglected. 
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6. SUMMARY 

The significance of our results for both primordial 3 He and chemical evolution of 3 He 
lies in the fact that the abundances we find are quite low. Hence we have paid particular 
attention to factors that might lead us to measure a lower abundance than that which is 
actually present. In this paper we conclude that it is unlikely that large ionization corrections 
are required for our high excitation state 3 He nebulae. We use the abundance to select 
for high excitation H II regions. We determined an ionization correction factor for the Paper 
I nebular sample. There would have to be a conspiracy between the density and ionization 
structure for these ionization corrections to be underestimated by a substantial amount. 
Our synthetic nebular models (§3]) do show that large ionization corrections are possible. 
Such models, however, require that there be rather specific density and ionization structure 
configurations. Our simple ionization correction fails when y\ changes significantly between 
the low density gas where 3 He + dominates the emission and the high density gas where 4 He + 
dominates. The more complex the source the larger our potential error in determining the 
ionization correction. We find both numerically and observationally that the simple sources 
give the best ionization corrections (Sj). 

For the sources that we classify as simple and which also have well determined line 
parameters, the ionization correction given in §H should be accurate. Indeed, the major 
uncertainty is likely to be the value adopted for y^GAL- Even for more complex sources our 
simple ionization correction should be reasonable for any nebula that does not have large 
fluctuations in excitation. We have shown this to be the case for the complex nebula W43 
where we have single-dish and interferometry data for both radio recombination line and 
continuum emission f@. 



For a subset of our 3 He sources, we use both 8 GHz and 18 GHz RRL data to derive 
y\. The two sets of observations probe different regions of the nebula, yet we in general find 
very similar values for y\. This is consistent with the notion that for our sources selected 
for high excitation there are not large ionization corrections. 

We can use this simple ionization correction analysis because unlike many of the light 
elements 3 He/H abundances accurate to ~ 10% can yield important conclu sions for both the 



primordial 3 He and the chemical evolution of 3 He ( IWilson fc Roodlll994j ). This is in stark 
contrast to the situation for 4 He where the 4 He/H abundance needs to be dete rmined to an 



accur acy of < 1% if one is to draw cosmologically significant inferences (e.g., Ilzotov et al. 



20071 ). Furthermore, because the majority of our 3 He nebulae are not optically visible it 
is not likely that we will ever have sufficiently detailed information about their physical 
properties to make a significally better correction. 
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Table 1. K-band (18 GHz) Radio Recombination Line Parmeters 
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Table 1 — Continued 
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Table 2. NEBULA Model Parameters a 
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a The distance to all model nebulae is 8 kpc. 
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Table 3. NEBULA Model Results 
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V = 8 GHz; HPBW = 3'.5. 
h v = 18 GHz; HPBW = 1(5. 
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Table 4. H II Region Excitation Properties 
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Table 4 — Continued 





Radio a 
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a Based on the Paper I X-band radio data that are appropriate for the 3 He + ionization 
correction. The 4 He + /H + abundance ratio is determined using the H91o; and He91a RRLs 
(only quality factors A and B are considered). The H-ionizing luminosity, N L , is calculated 
from the radio continuum emission. The spectral type is determined from A"l and a single 
star ZAMS model (Vacca et al. 1996; Panagia 1973, in parentheses). 

b Based on optical data from the literature. The references are listed in parentheses. 

References. — (1) Barlow et al. (1976); (2) Blum et al. (1999); (3) Chini & Wink (1984); 
(4) Conti & Alschuler (1971); (5) Conti & Leep (1974); (6) Crampton (1971); (7) Crampton 
et al. (1978); (8) Deharveng et al. (2000); (9) Esteban et al. (1999); (10) Feinstein & 
Vazquez (1989); (11) Fich & Silkey (1991); (12) Georgelin (1975); (13) Georgelin et al. 
(1973); (14) Goy (1980); (15) Haikala (1995) (16) Hanson & Conti (1995); (17) Hunter & 
Massey (1990); (18) Lynds & O'Neil Jr. (1986); (19) Miller (1968); (20) Moffat et al. (1979); 
(21) Moreno & Charvarria-K. (1986) (22) Peimbert et al. (1992b); (23) Persi et al. (1987); 
(24) Shaver et al. (1983); (25) Vflchez & Esteban (1996). 
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Fig. 1. — Spectra of the 91a X-band (8 GHz) and the 70a K-band (18 GHz) radio recombi- 
nation line emission for the H II regions NGC 7538 and W43. The vertical lines flag from left 
to right the C, He, and H transitions. (The additional flag in the 91a spectra is the H 154e 
transition.) The intensity scale is in units of milliKelvins. The LSR velocity increases from 
left to right at 1.3 km sec -1 per channel. 
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Fig. 2. — The 4 He + /H + , abundance ratio plotted as a function of the principal quantum 
number, n, of the RRL transition. The 70a, 88/3, 91a, and 114(3 transitions are shown. 
Abundances are plotted only for transitions that have a quality factor of A or B (see Table 
I and Paper I). 
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Fig. 3. — Comparison of the y\ ratios derived from the H91a and H70a spectra for the 10 
K-band sources in the order listed in Table [TJ Only abundances for quality factor A or B are 
plotted. Top Left: The H91a (circles) and H70a (diamonds) y^ ratios for each nebula. Top 
Right: The y^ (H91a) — yX (H70a) difference. Bottom Left: The fractional y\ deviation: 
(yt (H91a) -y\ (H70a)) j y\ (H91a). Bottom Right: The y\ (H91a) / yt (H70a) ratio. 
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Fig. 4. — The 4 He + /H + , y^, abundance ratio plotted as a function of the H-ionizing lumi- 
nosity, A^l- The yX values are calculated using only the H91a and He91a transitions. Only 
transitions that have a quality factor of A or B (see Paper I) are shown. Circles denote 
simple sources; triangles identify complex sources (see text). 



